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We study the various linear responses of neutron stars to external relativistic tidal fields. We focus 
on three different tidal responses, associated to three different tidal coefficients: (i) a gravito-electric- 
type coefficient Gjj.i — [length]^*"'"^ measuring the ^*'^-order mass multipolar moment GMa^...af^ 
induced in a star by an external ^"^-order gravito-electric tidal field Gai...ai \ (ii) a gravito-magnetic- 
type coefficient Gai = [length]^*^^ measuring the spin multipole moment GSai...ai induced in a 
star by an external £*'^-order gravito-magnetic tidal field Hai...a(\ and (iii) a dimensionless "shape" 
Love number hi measuring the distortion of the shape of the surface of a star by an external t^- 
order gravito-electric tidal field. All the dimensionless tidal coefficients G/if/i?^*^^, Goij B^^'^^ and 
hi (where R is the radius of the star) are found to have a strong sensitivity to the value of the 
star's "compactness" c = GM j (c^R) (where we indicate by co the speed of light). In particular, 
Gfj.t/ R^^^^ ~ ki is found to strongly decrease, as c increases, down to a zero value as c is formally 
extended to the "black-hole limit" (BH) c^^ = f/2. The shape Love number hi is also found to 
significantly decrease as c increases, though it does not vanish in the formal limit c —» c^^, but 
is rather found to agree with the recently determined shape Love numbers of black holes. The 
formal vanishing of fig and cr^ as c —> c^^ is a consequence of the no-hair properties of black holes. 
This vanishing suggests, but in no way proves, that the effective action describing the gravitational 
interactions of black holes may not need to be augmented by nonminimal worldline couplings. 

PACS numbers: 04.25.Nx, 04.40.Dg, 95.30.Sf, 



I. MOTIVATION AND INTRODUCTION 

Coalescing binary neutron stars are one of the most im- 
portant (and most secure) targets of the currently oper- 
ating network of ground-based detectors of gravitational- 
waves. A key scientific goal of the detection of the 
gravitational-wave signal emitted by coalescing binary 
neutron stars is to acquire some knowledge on the equa- 
tion of state (EOS) of neutron-star matter. Recent break- 
throughs in numerical relativity have given example of 
the sensitivity of the gravitational- wave signal to the EOS 
of the neutron stars [1-4]. However, this sensitivity is 
qualitatively striking only during and after the merger of 
the two neutron stars, i.e. for gravitational wave frequen- 
cies /gw ^ 1000 Hz, which are outside the most sensitive 
band of interferometric detectors. It is therefore impor- 
tant to study to what extent the gravitational- wave signal 
emitted within the most sensitive band of interferometric 
detectors (around /gw ~ 150 Hz) is quantitatively sensi- 
tive to the EOS of neutron stars. In such a regime, the 
two neutron stars are relatively far apart, and the prob- 
lem can be subdivided into three separate issues, namely: 

(i) to study the response of each neutron star to the tidal 
field generated by its companion; 

(ii) to incorporate the corresponding tidal effects within a 
theoretical framework able to describe the gravitational- 
wave signal emitted by inspiralling compact binaries; and 

(iii) to assess the measurability of the tidal effects within 
the signal seen by interferometric detectors. 

A first attack on these three issues has been recently 
undertaken by Flanagan and Hinderer [5, 6]. [See also [4] 
for an attempt at addressing the third issue.] Our aim 



in this work, and in subsequent ones, is to improve the 
treatment of Refs. [5, 6] on several accounts. The present 
work will focus on the first issue, (i), above, namely 
the study of the tidal response of a neutron star. Our 
treatment will complete the results of [6] in several direc- 
tions. First, we shall study not only the usually consid- 
ered "electric-type" , "tidal" , "quadrupolar" Love num- 
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ber Gii2 = -k2R'^, but also several of the other tidal 
o 

coefficients of a self-gravitating body. This includes not 
only the higher multipolar analogues G/x^ cx keR^^~^^ of 
/i2, but their "magnetic- type" analogues Gai ( first intro- 
duced in [7]), as well as their (electric) "shape- type" kin 
hg. Second, we shall study in detail the strong sensitivity 
of these tidal coefficients to the compactness parameter-'^ 
c = GM/cqR of the neutron star. Note, indeed, that the 
published version of Ref. [6] was marred by errors which 
invalidate the conclusions drawn there that fc2 has only a 
mild dependence on the compactness c (see e.g. Eq. (27) 
or Fig. 2 there). [These errors were later corrected in an 
erratum, which, however, did not correct Eq. (27), nor 
Fig. 2.] We shall interpret below the strong sensitivity of 
fif and cTi to c, and contrast the vanishing of /if and cri in 
the formal "black- hole limit" c 1/2, to the nonvanish- 
ing of the "shape" Love numbers hi in the same limit. In 
order to approach the "black-hole limit" (which is, how- 
ever, disconnected from the perfect-fluid star models), 
we shall particularly focus on the incompressible mod- 



^ To avoid confusion with the compactness, we sometimes denote 
the velocity of light as cq . 
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els which can reach the maximum compactness of fluid 
models, namely Cmax = 4/9. 

In subsequent works, we shall show how to incorpo- 
rate the knowledge acquired here on the various tidal 
responses of neutron stars into the Effective One Body 
(EOB) framework. Indeed, recent investigations [8, 9] 
have shown that the EOB formalism is the most accurate 
theoretical way of describing the motion and radiation of 
inspiralling compact binaries. 

This paper is organized as follows: Sec. II is an in- 
troduction to the various possible tidal responses of a 
neutron star. Section III discusses the relevant equations 
to deal with stationary perturbations of neutron stars 
that are then used in Sec. IV and Sec. V to compute 
the electric- type (fii) and magnetic- type (cr^) tidal co- 
efficients. Section VI is devoted to the computation of 
the "shape" Love numbers hi. Sections VII, VIII and IX 
provide explicit numerical results related to fit, ai and hi 
respectively. The concluding section. Sec. X, summarizes 
our main results. 



II. THE VARIOUS TIDAL RESPONSES OF A 
NEUTRON STAR 

Let us first recall that the motion and radiation of a 
system of well separated, strongly self-gravitating ("com- 
pact"), bodies can be theoretically investigated by a 
"matching" approach which consists in splitting the 
problem into two subproblems: 

(i) the outer problem where one solves field equations 
in which the bodies are "skeletonized" by worldlines en- 
dowed with some global characteristics (such as mass, 
spin or higher-multipole moments), and 

(ii) the inner problem where one obtains the near- 
worldline behavior of the outer solution from a study of 
the influence of the other bodies on the structure of the 
fields in an inner world tube around each body. 

This matching approach has been used: to obtain the 
dynamics of binary black holes at low post-Newtonian or- 
ders [10-12], to prove that the tidal deformation of com- 
pact bodies will start to introduce in the outer problem 
a dependence on the internal structure of the constituent 
bodies (measured by a "relativistic generalization of the 
second Love number" k) only at the fifth post-Newtonian 
(5PN) level [13], and to derive the dynamics of com- 
pact bodies in alternative theories of gravitation [14-16]. 
Finite-size corrections to the leading "skeletonized" dy- 
namics can be taken into account by adding nonminimal 
worldline coupHngs to the effective action [17, 18]. 

Let us start by considering the "inner problem" for 
a neutron star, i.e., the influence of the other bodies in 
the considered gravitationally interacting system ^. As 



explained, e.g. in Ref. [13], the matching method uses 
a multi-chart approach which combines the information 
contained in several expansions. One uses both a global 
weak- field expansion gfiu{x) = r]^,^ + h^i}}{x) + h^f^^(x) + . . . 
for the outer problem, and several local expansions of the 
type 



gZ in) 



[xD + 



(1) 



for each inner problem. Here, G^^^ denotes the metric 
generated by an isolated neutron star, as seen in a local 
inner coordinate system X^, which is nonlinear ly related 
to the global ( "barycentric" ) coordinate system by an 
expansion of the form 



1 



(2) 



Here, the suffix A = 1,2,...,N labels the considered 
member of the TV-body system, while H^^ [X]^) denotes 
the metric perturbation, seen in the local A-frame, Ba- 
cause of the combined influence of the various compan- 
ions B ^ A oi A. In the leading approximation is a 
sum of separate contributions due to each B ^ A: each 
contribution then contains both the far-away field gener- 
ated by the B worldline, its deformation as it propagates 
on the "background" metric c'fl generated by A, and 
the tidally-induced effect of the deformation of A by the 
effect of B. 

Before tackling the technical problem of computing 
H^p, let us recall the general structure of tidal expan- 
sions in general relativity [7, 12, 19]. We will use here the 
notation and results of the general multi-chart approach 
to the general relativistic dynamics of N self-gravitating, 
deformable bodies developed by Damour, Soffel and Xu 
(DSX) [7, 19-21]. 

Using the DSX notation (with T = X^jc), 



Go1,(X) = -exp(-2M^^/c2) 



r3 a J 



E^{X) = daW 



B^{X) = eabcdb {-m^ 



(3) 
(4) 

(5) 
(6) 



one defines, in the local frame of each body A, two sets 
of "gravito-electric" and "gravito-magnetic" relativistic 
tidal moments, and i/^, respectively as^ 



Ht{T)^\d^L-iBl)iT.y^) 



(7) 
(8) 



^ In the following, we shall have in mind a binary system made 
cither of two neutron stars or of a neutron star and a black hole. 



^ As in DSX, L denotes a a multi- index a\a2 ■ ■ ■ ai and (ai . . .ai) 
a symmctric-tracc-frec (symmetric-trace-free) projection. 
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where and denote the externally-generated parts 
of the local gravito-electric and gravito-magnetic fields 
and B^. In the presently considered approxima- 
tion where G'^j^'^ is stationary, and where it is enough 
to consider the linearized, multipole expanded, perturba- 
tion H^^p in Eq. (1), the externally generated parts E^ 
and B^ are well-defined and capture the terms in E^ and 
B^ that asymptotically grow as R^~^ as i? = |X| oo. 
The (seemingly contradictory) formal limit X° in- 
dicated in Eqs. (7)- (8) refers to the matching performed 
in the outer problem (where, roughly speaking, the outer 
limit J^^outcr ^ can still refer to a worldtube which is 
large, in internal units, compared to the radius of body 
A). 

Besides the externally-generated "tidal moments" (7)- 
(8), one also defines the internally-generated "multi- 
pole moments" of body A, M^{T) (mass moments) and 
S^{T) (spin moments) as the symmetric-trace- free ten- 
sors that parametrize the body-generated terms in the 
metric coefficients W^, that asymptotically decrease 
(in the ^-body zone) as as i? = |X| ^ oo. The 

normalization of these quantities is defined by Eqs. (6.9) 
of [f 9] (and agrees with the usual one in post-Newtonian 
theory) . 

In the stationary case (which is relevant to our present 
"adiabatic" approach to tidal effects), this normalization 
is such that the "internally-generated" post-Newtonian 
metric potentials W+^{X), W+'^iX) read 



(9) 



e>i 



R 



-^da (A^ - A) , 



(10) 



where A'^ — A is a gauge transformation (which would 
drop out if we had considered the gravito-magnetic field 

In the first post-Newtonian approximation considered 
by DSX, the separation of the local-frame potential 
W^{X) into an "externally-generated" part and an 
"internally-generated" one W^^, is well defined (thanks 
to the structure of Einstein's equations). In the case 
we are considering here of a linearly perturbed, quasi- 
stationary, fully relativistic neutron star, the asymptot- 
ically growing character (as i? — > oo) of the externally- 
generated potentials allows one to uniquely define the 
tidal moments (7)-(8). On the other hand, the asymp- 
totic decrease cx of the internally generated 
multipolar potentials (9)-(10) introduces an ambiguity 
in their definition. For an attempt to uniquely de- 
fine the gravito-electric quadrupole moment in- 
duced on a black hole by an external tidal moment 
see [22]. Here, instead of relying on such a conventional 
(harmonic-coordinates related) definition of the induced 



multipole moments, we shall follow the spirit of Sec. 5 of 
[13] in defining Mf, Sf^ as parametrizing the (uniquely 
defined) pieces in the local- frame metric G^p{X) which 
violate the "effacing principle" , in that they directly de- 
pend on the body A being a neutron star, rather than a 
black hole. Reference [13] explicitly treated the dominant 
even-parity case, and introduced (see Eq. (11) there) a 
"dimensionless constant fc" (= 02 as defined below) as a 
"relativistic generalization of the second Love number". 
This minimal definition (which will be made fully precise 
below) is rather natural, and coincides with the definition 
adopted in [5, 6]. 

With this notation in hands, we can define the two 
"tidal-polarizability" coefficients /i£ and at introduced in 
Eqs. (6.19) of [7]. These coefficients relate the (electric 
or magnetic) tidal induced multipole moments to the 
corresponding external tidal moments, i.e. 



Mi 



(11) 
(12) 



The electric-type (or "even-parity") tidal coefficient fif 
generalizes the (fc^-type) Newtonian "Love number" . For 
the leading quadrupolar tide, /i2, as defined by Eq. (11), 
agrees with the quantity denoted A in [5, 6]. The 
magnetic-type (or "odd-parity" ) quadrupolar tidal co- 
efficient (T2 is proportional to the quantity 7 which has 
been considered in the investigations of Favata [23] which 
were, however, limited to the first post-Newtonian ap- 
proximation. Here, we shall consider the case of strongly 
self- gravitating bodies (neutron stars), and study the 
dependence of both and af on the compactness 
CA = {GM/c^R)a of the considered neutron star. Let 
us also note that, in terms of finite-size corrections to 
the leading point-particle effective action Spointmass = 
— Y^a! ^AdsA, the two tidal effects parametrized by 
fi( and ai correspond to nonminimal worldline couplings 
respectively proportional to 



J dsA {Gif , and 4 J dsA {HfY 



(13) 



The leading, quadrupolar corrections (13) can be repro- 
duced (using the link between G^j, and u'^UA^^aub^ ^^^^ 



Hab and 



-u'XR^acd^ see Sec. 3.D of [19] and [12]) as 
the following nonminimal couplings involving the Weyl 
tensor 



J dsEc^pE"'^ and G2 j dsBo^pB'"^ 



(14) 



where 
sors £, 



af3 



dz'^/ds, and we have introduced the ten- 
= u''u''Gfj,a„l3 and 8^/3 = Wn'^C* a, with 



* Here, we consider a nonrotating star which is spherically sym- 
metric (with vanishing multipole moments) when it is isolated, 
so that and represent the multipole moments induced 

by the influence of the external tidal fields and H^. 
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C^!/a/3 = h'^ P'yap being the dual of the Weyl ten- 
sor. In = 3 -|- 1 dimensions, and in absence of parity- 
violating couplings, the two terms (14) arc the only pos- 
sible isotropic couplings. In higher dimensions, there are 
three nonminimal isotropic couplings quadratic in the 
Weyl tensor as indicated in Eq.(90) of [18]. Note that we 
are using here the freedom of locally redefining the dy- 
namical variables to eliminate terms proportional to the 
(zeroth-order) equations of motion, such as terms involv- 
ing the Ricci tensor; see, e.g., the discussion of finite-size 
effects in tensor-scalar gravity in Appendix A of Ref. [17]. 

Let us finally note that there are other "tidal coef- 
ficients" which might be interesting to discuss. First, 
though the linear relations (11)-(12) arc the most gen- 
eral ones that can exist in the (parity-preserving) case 
of a nonspinning neutron star, the tidal properties of 
a spinning neutron star will involve other tidal coeffi- 
cients, proportional to the spin, and associated to a mix- 
ing between electric and magnetic effects. Such electric- 
magnetic mixing terms would correspond, say in the lead- 
ing quadrupolar case, to nonminimal worldline couplings 
quadratic in C^j/q^ and linear in the spin tensor S"^jj. 

There exist also other tidal coefficients which do not 
have a direct dynamical meaning, but which generalize 
the "first type" of Love numbers introduced in the theory 
of Newtonian tides. Indeed, it is physically meaningful to 
define, for any a "shape" Love number measuring the 
proportionality between the external tidal influence, and 
the deformation of the geometry of the surface of the con- 
sidered (neutron) star. More precisely, limiting ourselves 
to the electric-type tides, one can define a dimcnsionless 
number by writing, as one does in Newtonian theory. 



disturb 



(15) 



or, equivalently. 



5R 
~R 



disturb 



gR 



(16) 



where {6R/R)i oc Pe{cos9) represents the fractional de- 
formation of the (areal) radius R of the neutron star 
(measured in a geometrically invariant way, by relat- 
ing it to the inner geometry of the deformed surface), 
where J7^^"'*"'''^(i?) cx R^Pg{cos9) represents the usual, 
external, Newtonian tidal potential deforming the star, 
formally evaluated at the radius of the star (as if one 
were in fiat space), and where g = GM/R^ represents 
the usual Newtonian surface gravity of the neutron star. 
This hi, "shape" Love number has been recently consid- 
ered in the theory of the gravitational polarizability of 
black holes [24] and it will be interesting to compare and 
contrast the values of the hi for black holes to the values 
of hi for neutron stars, especially in the limit where the 
compactness gets large. See Section VI below which will 
give the exact definition of the quantity (SR/R)^. 



III. STATIONARY PERTURBATIONS OF A 
NEUTRON STAR 



The unperturbed structure of an isolated (nonrotating) 
neutron star is described by a metric of the form 

-e-^Mdi^ + e^'^^'Ur^ + r^d9?. (17) 



Here, and in the following, for notational simplicity we 
shall denote the local (spherical) coordinates of the A- 
body frame simply as {t, r, 6, ip) (with dft^ = d6^ + 
sin^ 6d(p^), instead of the upper case letters (T, R, 6, $) 
that would more closely follow the DSX notation recalled 
above. Introducing as usual the radial dependent mass 
parameter m(r) by^ 



1 



2m(r) 



(18) 



and assuming a perfect-fluid energy-momentum tensor 

= {e + p)u^Ui, +pg^,u, (19) 

the spherically symmetric metric coefficients v{r), m{r) 
and the pressure p{r) satisfy the Tolman-Oppenheimer- 
Volkoff (TOV) equations of stellar equilibrium 



dm 2 
— — = 47rr e, 
dr 



dp 
dr 



, , m -|- ^-KT^p 



dv 1(m, + 47rr'^p) 
dr r^ — 2mr 



(20) 
(21) 
(22) 



These equations are integrated from the center outward 
once that a barotropic EOS relating p to e is pro- 
vided. We shall consider several types of barotropic 
EOS, namely two different types of "polytropic EOS" 
(an e-polytrope, with p = kc'' , and a /i-polytrope" , with 
p — kfi^ and e = fi + p/(7 — 1), where fi — nmi, is the 
baryonic rest-mass density), and two different tabulated 
("realistic") EOS (the EPS one [25] and the SLy [26] one). 
In view of the current large uncertainty in the correct 
description of dense nuclear matter, we are not claiming 
that our selection of "realistic" EOS is physically pre- 
ferred (see, e.g., Ref. [27] and references therein for a 
thorough comparison among models from various EOS). 
We have chosen them because they have been used in 
recent numerical relativity simulations of binary neutron 
star coalescence [28, 29]. As for the polytropic EOS, they 
have also been often used in numerical relativity simula- 
tions (especially the ^-polytrope one), and their depen- 
dence on the adiabatic index^ is a convenient way of vary- 
ing the "stiffness" of the EOS (the limit 7 ^ oo repre- 
senting the stiffest possible EOS, namely incompressible 
matter with e = const and an infinite speed of sound). 



Henceforth, we shall often set G = c = 1 
® As is well-known, the dependence on the "polytropic constant" k 
can be absorbed in the definition of suitable "polytropic units" . 
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BacausG of the spherical symmetry of the background, 
the metric perturbation 

G^piX) = Gt^°\x)+H^0iX), (23) 

here considered at the hnearized level, can be analyzed 
in (tensor) spherical harmonics. The metric is expanded 
in even-parity and odd-parity tensor harmonics as 

(24) 

In the Reggc- Wheeler gauge, and following standard def- 
initions for the expansion coefficients and the sign con- 
ventions of [30; 31], one has 

+ H^'^e^dr^ + r^K'^"'dn^] Y«™, (25) 

while the nonvanishing components of H^^^ are H^'^ = 

hoe/VeYim and H^°2 = hie/VsYi^ where {A,B) = 
(6*, tp) and where is the mixed form of the volume 
form on the sphere S^. 

Our aim is then to solve the coupled system of the per- 
turbed Einstein's equations, together with the perturbed 
hydrodynamical equations V"^TQ^[e,p] = 0, so as to de- 
scribe a star deformed by an external tidal field. We 
shall only consider stationary perturbations ("adiabatic 
tides"). 

A. Even-parity, stationary barotropic 
perturbations 

Even-parity, stationary perturbations of a barotropic 
star simplify in that: (i) the metric perturbations reduce 



to two functions H = Hq = H2, and K (with Hi = 0), 

(ii) the fluid perturbations are described by the logarith- 
mic enthalpy function h, such that Sh = 6p/{e +p), and 

(iii) the latter logarithmic enthalpy function is simply re- 
lated (in absence of entropy perturbation) to the metric 
function H by 



6h = H. 

2 



It was then showed by Lindblom, Mendell and Ipser [31] 
how to convert the system of first-order radial differential 
equations relating H' , K', H and to a single second- 
order radial differential equation for the metric variable 
H (such that Hqq = —e" HYim) of the form 



H" + CiH' + CqH = 0. (27) 



[As usual, we shall generally drop the multipolar index £ 
on the various metric perturbations. The presence of a 
factor Yirn{0,'p), or Pf(cos6'), in (or to be added to) the 
considered metric perturbation is also often left implicit.] 
Taking the stationary limit [uj —^ 0) of the results given in 
Appendix A of [31] (together with the barotropic relation 
SU ~ 0) one gets 



Ci = - + i(z.'-A') 
r 2 



— + 47rr(p — e) 



Co = 



i(e+l) , ,de 
- ' +47r(e+p — -f 47r e+p) 

dp 

((i+l) . . .de 



+ v" + [y'f + (2 - rv') {?>v' + A') 
Ir 



+ 47r(e + p) — + 47r (5e + 9p) 
dp 



I '\2 



(28) 



(29) 



where we have used the background (TOV) equations to 
rewrite Ci and Cq. As a check, we have also derived from 
scratch Eq. (27) by starting from the "gauge-invariant" 
formalism of Ref. [32] . Equation (27) generalizes to an 
arbitrary value of the multipolar order £ Eq. (15) of 
Ref. [6], which concerned the leading quadrupolar even- 
parity tide. 

For completeness, let us note that the other metric 
variable, can be expressed as a linear combination of 



H and H' , namely 

K = aiH' + a2H, (30) 

where the explicit expressions of the coefficients ai and 
0:2 can also be deduced by taking the stationary limit of 
the results given in Appendix A of [31]. 
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B. Odd-parity, stationary perturbations 

It was shown by Thorne and Campolattaro [33] that 
odd-parity perturbations of a nonrotating perfect-fiuid 
star consists only of metric fluctuations, and do not affect 
the star's energy density and pressure. One might naively 
think that this means that an odd-parity tidal field will 
induce no (gauge-invariant) spin multipole moment in a 
(nonrotating) star. This conclusion is, however, incorrect 
because the "gravitational potential well" generated by 



the stress-energy tensor of the star does affect the "radial 
propagation" of the external odd-parity tidal fields and 
necessarily adds an asymptotically decreasing "induced" 
tidal response to the "incoming" tidal field. To describe 
this phenomenon, it is convenient to describe the odd- 
parity perturbation by means of the (static limit of the) 
"master equation" derived by Cunningham, Price and 
Moncricf [34] (see also Rcf. [35]). In the stationary limit, 
and in terms of the ordinary radial variable r (rather than 
the "tortoise" coordinate r,) this equation reads 



J 



i)" + ^[2m + Anr^ip - e)] ip' - 



1) 



6™ , / X 
— + 4^(e - p) 



xjj = 0. 



(31) 



In terms of the variables {ho, hi) entering the odd-parity 
perturbations, the odd-parity master function can 
be taken to be either e'^'^~'^^/'^hi/r, or the combination 
rdthi — r^driho/r^) [see e.g. [36] for more details]. As hi 
vanishes in the stationary limit, we can define ^ as being 



= r^dr 



rhg — 2/7,0. 



(32) 



IV. 



COMPUTATION OF THE ELECTRIC-TYPE 
TIDAL COEFFICIENT m 



The electric- type tidal response coefficient fii, defined 
by Eq. (11) above, can be obtained by going through 
three steps: (i) numerically solving the even-parity mas- 
ter equation (27) within the neutron star; (ii) analytically 
solving the same master equation (27) in the exterior of 
the star; and (iii) matching the interior and exterior solu- 
tions across the star surface, taking into account the defi- 
nition (11) to normalize the ratio between the "growing" 
and "decreasing" parts of H{r), namely 7yg™™'is ~ 
versus - /x^r-(^+i). 



A. The internal problem 

The internal value of the metric function H{r) is ob- 
tained by numerically integrating Eq. (27) , together with 
the TOV equations (20)- (22), from the center (or, rather 
some very small cut-off radius rg = 10~^) outwards, 
starting with some central values of to,p, i^, H and H' . 
For iJ, one takes as starting values at the cut-off radius 



H{ro) = and i/'(ro) 



The latter boundary 



conditions follow from the analysis of Eq. (27) around 
the regular-singular point r = 0, which shows that 
H{r) ~ hr^ (where h is an arbitrary constant) is the 
most general regular solution around r = 0. As Eq. (27) 
is homogeneous in H , the scaling constant h is irrelevant 



and will drop out when we shall match the logarithmic 
derivative 



y 



rH' 



(33) 



across the star surface. This is why it is enough to use 
h = 1 as initial boundary conditions for H . 

The main output of this internal integration procedure 
is to compute (for each value of €) the value of the internal 
logarithmic derivative (33) at the star's surface, say r = 
R 



B. The external problem 



(34) 



As noticed long ago by Regge and Wheeler [37] and 
Zcrilli [38], the exterior form of the stationary, even- 
parity master equation (27) (e = p = 0, m(r) = M) 
can be recast as an associated Lcgendre equation (with 
1 = 1 and m = 2). More precisely, in terms of the inde- 
pendent variable x = r/Ad — 1, the exterior form of (27) 
reads 



{x^ - l)H" + 2xH' - £{1 + 1 



X — 1 



H = 0, (35) 



where the prime stands now for d/dx. Its general solution 
can be written as 



H = apPi>2{x) 



n 



(36) 



where the hat indicates that the associated Legendre 
bmctions of first, Pi2, and second^, kind have been 



Note that, contrary to the usual mathematical definition of 
Qimi^), which is tuned to the real interval — 1 < a; < we 
need to work with Qij^{x) in the interval a; > 1, This means 

replacing log ( — — — j with log ( ~^ 



1 



1 
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normalized so that Qi2 — ~ [M/rY'^^ and 

Pi2 ~ ~ {t/MY when x ^ oo or r ^ oo\ aq and 
ap arc integration constants to be determined by match- 
ing to the internal solution. Defining = ag/ap, the 
exterior logarithmic derivative y°^^ = rH' / H reads 



yr\x) 



{l + x) 



Pg2{x) + aiiQ'i,2{x) 
Pt2{x) + atQnix) 



(37) 



Matching at the star's surface, and computation 
of the "electric" tidal Love number 



As Eq. (27) is second-order in the radial derivative of 
H , one expects that H and H' will be continuous at the 
star's surface. Actually, the issue of regularity at the star 
surface is somewhat subtle because some of the thermo- 
dynamic variables (such as pressure) do not admit regular 
Taylor expansions in r — i? as r ^ i?. For instance, while 
the logarithmic enthalpy h{p) — j^dp/{e-]-p) vanishes 
smoothly {h{r) oc r—R) across the surface, one finds that 
(for any polytrope) p[r) oc (r — i?)T/(T~i) and that the 
term involving the inverse of the squared sound velocity 

= dp/de in Eq. (27) is singular (when 7 > 2). namely 



dp 



rY 



(38) 



Despite this mildly singular behavior of the coefhcicnt 
Co of (27) and despite the fact that the exact location 
of the tidally-deformed star surface is slightly displaced 
from the "background" value r = R, one checks that it 
is correct (when 7 < 00) to impose the continuity of H 
and H' aX r = R. [Note that we consider here the case 
of a finite adiabatic index 7. The incompressible limit 
7 — > 00 leads to a master equation which is singular at 
the surface, and which must be considered with care. See 
below our discussion of the incompressible limit.] This 
continuity then imposes the continuity of the logarithmic 
derivative rH' /H . This leads to the condition y°^^{R) = 
?/™*(i?) = yg, which determines the value of the ratio 
ag ~ aq/ap in terms of the compactness c = M / R of the 
star 



ae. 



P'(,2{x) - cyiPi2{x) 



(39) 



On the other hand, the ratio = aq/ap can be related 
to the tidal coefficient p,i by comparing (modulo an over- 
all factor —2), 



(SHooe 



growing 



^growing ^^-j 



(40) 



= apPi2{x)Yi„i c:i ap l^—j Yim, 

-z^j decreasing ^decreasing 

/ r 

= aqQ£2{x)Yim ~ aq Yi„i, (41) 



respectively to 

W 
W+ 



{-Y 



1 



L I 



(42) 
(43) 



(see e.g., Eq. (4.15a) of Ref. [7]) where n" = X°- /r is a 
radial unit vector. Using the fact that 

dLT-^ = (-)^(2£- l)!!ri^r~(^+i\ (44) 

and M]^ = fifGL, and remembering that Glu^ cx 
y«m(6', '^)^ we see that 



i2i- 



(45) 



Note that Gpi has the dimensions of [length] There 
are then two natural ways of expressing G/i^ in terms 
of a dimensionless quantity. Either by scaling it by the 
(2£ + l)-th power of GM/cl, which leads to 



(2£- 1)!! 



G^JL|^ 



{GM/4 



2\2e+i 



ae, 



(46) 



or by scaling it by the (2£+ l)-th power of the star radius 
R, which gives 



(2^-1)!! 



Gfii 



■^21+1 

Alternatively, we can write 



2ki — aiC 



2e+i 



Gfii 



ai /GM\ 



2i+\ 



2ko 



(2^-1)!! 



R 



2£+l 



(47) 



(48) 



The scaling of G/x^ by means of R?'^^^ is the traditional 
"Newtonian" way of proceeding, and leads to the intro- 
duction of the dimensionless "second tidal Love number" 
kg (conventionally normalized as in Eq. (47) above). 
One can finally write kg as 



k, = 4c2^+^ 



1 2i+\ Pki^) - cyiPn{x) 
ai ~ c — 

2 Q'i2{x) - cyiQnix) 



x = l/c-l 

(49) 

The dimensionless Love number kg has the advantage 
of having a weaker sensitivity on the compactness c = 
GM/{cqR) (especially as the compactness formally tends 
to zero, i.e. in the Newtonian limit). Note, however, that 
the dimensionless quantity which will most directly enter 
the gravitational-wave phase of inspiralling binary neu- 
tron stars (NS) is Gne/{GM/c^oY^+^ ^ ag ^ c-'^^^+^^kg. 

The evaluation of the result (49) for kg yields the fol- 
lowing explicit expressions for 2 < ^ < 4 (with, for sim- 
plicity, y = ye): 
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fc2 = 7(l-2c)V[2c(y-l)-y + 2] 
5 



|2c (4(?/ + l)c^ + (6y - 4)c3 + (26 - 22yy + 3(5y - 8)c - 3?/ + 6) 
3(l-2c)2(2c(2/-l)-z^ + 2)log 



1 - 2c 



fca = -(1 - 2c)2c^ [2{y - l)c2 - 3(2; - 2)c + y - 3] 



(50) 



X |2c [4(2/ + l)c^ + 2(9y - 2)c'' - 20(7?/ - 9)c^ + 5(37y - 72)c2 - 45(2y - 5)c + 15{y - 3)] 
- 15(1 - 2c)2 (2(y - l)c2 - 3(2/ - 2)c + y - 3) log ( } ' 



(51) 



32 

A;4 = Yi^^l - 2c)2c9 [I2(y - l)c3 - 34(2/ - 2)c2 + 28(2/ - 3)c ~ 7{y ~ 4) 



X <{ 2c [8(2/ + l)c^ + (682/ - 8)c5 + (1284 - 9962/)c'' + 40(552/ - 116)c^ + (5360 - I9l0y)c^ + 105(7?/ - 24)c - 105(2/ - 



15(1 - 2c)2 [12(2/ - l)c3 - 34(2/ - 2)c2 + 28(2/ " 3)c - 7(2/ - 4)] log J^- 

I 



1 



Equation (50) above agrees with the corrected version of 
Eq. (23) of [6] . Note that, independently of the values of y 
(as long as it does not introduce a pole singularity, which 
will be the case), the results (50)-(52) (and, more gener- 
ally, the result (49)) contain an overall factor (1 — 2c) ^ 
which formally tends (quadratically) to zero when the 
compactness c = GM/R "tends" toward the compact- 
ness of a black hole, namely c^" = GM/{2GM) = 1/2. 
[The singular logarithm log[l/(l — 2c)] in the denomi- 
nator is also easily checked to be always multiplied by 
(1 — 2c) ^ and thereby not to affect the cx (1 — 2c) ^ for- 
mal vanishing of kg as c 1/2.] This property can 
be easily understood as a consequence of the "no-hair" 
properties of black holes. Indeed, among the two so- 
lutions of the exterior tidal perturbation equation (27), 
the no-hair property means that the solution which is 
"rooted" within the horizon, i.e., the "asymptotically de- 
creasing" solution Qe2{x) is singular at the horizon, i.e. 
when X = R/M — 1 ^ 1. More precisely, this singular 
behavior is 



Ql2{x) 



/ 2 _ ,,2/2 d'^Qt{x) 



(x' 



dx^ 



log 



X - 1 



Peix) 



(x^ ^ l)ix - 1)-^ ^ {x - 



(53) 



so that the most singular term in the denominator of ae 
or h is Q'e2{x) - {x - 1)^^ _ (ij _ 2M)~^ - (1 - 2c)-^ 
which is at the origin of the presence of a factor (1 — 2c)^ 
in a£ and fc^. One might naively think that this behav- 
ior proves that the "correct" value of the k£ tidal Love 
numbers of a black hole is simply kf^ = 0. However, 



we do not think that this conclusion is warranted. In- 
deed, as we explained above, the definition used here (and 
in [5, 6, 13]) of the Love numbers of a (neutron) star 
consists in selecting, within the gravitational field of a 
tidally distorted star, the terms which violate the "effac- 
ing principle" (in the sense of Ref. [13]), i.e. the internal- 
structure-dependent terms which differentiate the tidal 
response of a (compact) star, from that of a black hole. 
From this point of view, the vanishing of fc^ as c ^ c^^ 
is mainly a consistency check on this formal definition. 
The question of computing the "correct" value of kg for 
a black hole is a technically much harder issue which 
involves investigating in detail the many divergent dia- 
grams that enter the computation of interacting point 
masses at the 5-loop (or 5PN) level. 

Indeed, the issue at stake is the following. When de- 
scribing the motion of two black holes (as seen in the 
"outer problem") by a skeletonized action of the form 

S = 5'pointmass + S'nonminimal, thc prCSCUCe of UOnmiu- 

imal worldline couplings S'nonminimal of the type (13) 
and (14) can only be detected if one treats (when us- 
ing perturbative expansions in powers of G) the gen- 
eral relativistic nonlinear self-interactions entailed by 
Spointmass = - J2a I ^ Ads A at the ordcr of approxi- 
mation corresponding to S'nonminimal- For a black hole 
(of "radius" Ra = 2GMa/cq), the leading nonmini- 
mal coupling parameter scales as /i^ ^ fc^i?^/G ^ 

Rapjs « GMb) the 



k^G'^M\, so that (using 

leading nonminimal interaction /x^ J dsA^^pS^" is pro- 
portional to fc^G^M^Mg. The presence of an overall 
factor G^ (which is the same factor G^ that appeared 
in Eq. (19) in Sec. 5 of [13]) signals that such an ef- 
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feet is smaller than the leading (Newtonian) inter- 
action (cx GMaMb) between two point masses, so that 
it corresponds to the 5PN level. In the diagrammatic 
language of (post-Minkowskian or post-Newtonian) per- 
turbation theory (as used, e.g., in [39]), this corresponds 
to the 5-loop level. Let us recall that the computation 
of the interaction of two black holes at the 3-loop level 
was a technically complex enterprise that necessitated 
the careful consideration of many divergent diagrams, 
and the use of the efficient method of dimensional regu- 
larization [40, 41]. At the 3-loop level the result of the 
computation was (essentially) finite, though the use of 
harmonic coordinates in one of the computations [41] 
introduced some gauge-dependent infinities. As argued 
long ago [13], and confirmed by an effective action ap- 
proach [18], one expects to sec real, gauge- independent 
infinities arising at 5-loop (5PN), i.e. at the level where 
the effacing principle breaks down, and where, as ex- 
plained above, a parameter (~ ^2) linked to the internal 
structure of the considered compact body starts to en- 
ter the dynamics. Until a careful analysis of the 5PN 
nonlinear self-interactions is performed, one cannot con- 
clude from the above result (fc^^ — > as c ^ c^^) that 
the effective action describing the dynamics of interact- 
ing black holes is described by the pure point-mass action 
— J TTiAdsA without the need of additional nonmini- 
mal couplings of the type of Eq. (14). 

We have phrased here the problem within stan- 
dard (post-Minkowskian or post-Newtonian) perturba- 
tion theory, because this is the clearest framework within 
which the issue of higher order nonlinear gravitational 
interactions of point masses is technically well defined 
(when using, say, dimensional regularization to define the 
pcrturbative interaction of point masses in general rela- 
tivity [40; 41]). Note that, in the extreme mass ratio limit 
{Ma ^ Mb), where one might use black hole perturba- 
tion theory, the interaction associated to the leading non- 
minimal coupling parameter /i^ of Ma is proportional to 
M^ (see above). This is well beyond the currently stud- 
ied "gravitational self-force" effects, which are propor- 
tional to M\, and correspond to a "1-loop" effect within 
a black hole background. 



V. COMPUTATION OF THE MAGNETIC-TYPE 
TIDAL COEFFICIENT at 

The magnetic- type tidal response coefficient ai, de- 
fined by Eq. (12) above, can be obtained by following 
three steps, which are similar to those followed for the 
electric- type coefficient fii. 



now obtained from the behavior tp oc r^~^^ of the general 
regular solution at the origin. Again, the main output of 
the internal integration procedure is to compute (for each 
value of £) the value of the internal logarithmic derivative 
of ip, at the star surface, say 



V'ii 



(54) 



r=R 



B. The external problem 

As noticed long ago by Regge and Wheeler [37], the 
stationary odd-parity perturbations can be analytically 
solved in the exterior region. Similar to the even-parity 
case there exist two types of exterior solutions: a "grow- 
ing" type solution, say ipp{f), with f = r/M, and a "de- 
creasing" type one, say ipQi^)- We normalize them so 
that '4'p{t) — f^^^, and '!/'q(^) — as f ^ 00. The 
general analytical forms of t/jp and ij^q, for any can 
be obtained from Ref. [37]. In the case of the leading 
quadrupolar odd-parity perturbation, i = 2, the "grow- 
ing" analytical exterior solution of (31) is the very simple 
polynomial 



(55) 



while the "decreasing" one can be expressed in terms of 
an hypergeometric function F{a, b; c; z) as 



F 1,4;6; 



(56) 



The normalization of TpQ{f) is such that ipgif) ~ as 
r ^ 00. Note also that, for the special values a = 1, 
6 = 4, c = 6, the hypergeometric function is actually 
expressible in terms of elementary functions. The result 
has the form 

i/'Q=^(f) = Asf^ log (^!—-^^+A2f'^+Aif + An + A^if~\ 

(57) 

Going back to an arbitrary i, the general exterior solution 
of Eq. (31) can be written, analogously to the even-parity 
Eq. (36), as 



^-t = 6pV'p(r) + 6Q^Q(f). 



(58) 



This result allows one to compute the logarithmic deriva- 



tive y' 



,odd 



rip' /xj} of ijj in the exterior domain, namely 



yodd(^) = r 



>i.(f) + &,^' (f) 



ipp{f) + biipQ{f) ' 



(59) 



where be = bg/bp. 



A. The internal problem 

The internal value of the odd-parity master function tp 
is obtained by numerically integrating Eq. (31), together 
with the TOV equations. The boundary conditions are 



C. Matching at the star surface, and computation 
of the "magnetic" tidal Love number 

We again impose the continuity of ip, ip' , and therefore 
yodd _ j.^' j^^ at the star's surface. Similarly to the 
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even parity case, this determines the value of the ratio 
bi = bq/bp in terms of the compactness of the star: 



bf^- 



(60) 



f— 1/ c 



Again, we see at work the effect of the "no-hair" property 
in that the term ^'qif) in the denominator of (60) will 
become (from (57)) singular as (f — 2)~^ when f — > 2. 
This implies that bi will vanish proportionally to the first 
power of 1 — 2c in the formal limit where the star's com- 
pactness c c^^ = 1/2. 

The dimensionless quantity 6^, Eq. (60), is the odd- 
parity analog of the even-parity quantity ai, Eq. (39). 
In the even-parity case, was, essentially, the tidal re- 
sponse coefficient G^n scaled by {GM/cqY^^^ . In the 
present odd-parity case, the tidal response coefficient Gat 
has again the dimension [length]^^+^, and bi (for a gen- 
eral £) is essentially bi - G(Jt{GM / cl)-^'^^+^\ Before 
working out the exact numerical coefficient in this pro- 
portionality, we can note that the odd-parity analog of 
ki (i.e. essentially ki ~ G^i/R'^^^^ ^ c^^^^ai) will be 
obtained by scahng Gae by the {2i + l)-th power of the 
star radius R, and will therefore involve the new dimen- 
sionless combination 



Ji 



tp'p{r) - cyoddfppif) 



V'q(?^) - cyoddil'Qi'r) 



(61) 



--1/C 



One expects that the new odd-parity dimensionless com- 
bination ji (61) will, like ki, depend less strongly on the 
value of the compactness c than bi itself. 

Let us now derive the precise link between the odd- 
parity tidal response coefficient Gcr^, defined by Eq. (12), 
and the dimensionless quantities bt, Eq. (60), or 6fC^^+^, 
Eq. (61). To relate them, we start by noting that 
the Regge- Wheeler metric function /iq entering the odd- 
parity master quantity -0, Eq. (32), parametrizes the 
time X angle off-diagonal component of the metric per- 
turbation 



(62) 



where A, i? = 2, 3 = 0, iy9 are indices on the background 
coordinate sphere of radius r. The metric on S*^ 
is jAsdx^dx^ = r^dfl^, while = 7^*^eAc denotes 
the mixed form of the volume form ^eAsdx^ A dx^ ~ 
sin 9d9 A dip on S^. Let us now consider the gravito- 
magnetic field Ba, as defined by DSX. Modulo an ir- 
relevant numerical factor, it is the 3-dimensional curl 
of the time-space off-diagonal metric component: Ba oc 
^abcdbHoc- Let us focus on the "radial component" of the 
gravito- magnetic field Ba, i.e. the pseudo-scalar 



n • B = n^Ba (x n^EabcdbHa, oc e^^^AH, 



OB 



(63) 



Using then Eq. (62), one finds that 
n • B cx -ho{r)e'^^eJ^\/B^cYe„, = 



imr,. 



(64) 



where one used e^^ejf — j^'^' and the fact that 
■^^^S/A^BYem = + l)r-^Yi,n, where the factor 

comes from the fact that jab is the metric on a 
sphere of radius r, rather than a unit sphere. [The vari- 
ous proportionality signs refer to irrelevant, coordinate- 
independent, numerical factors.] Finally, we have the 
link 



^ = r^dr (^) cxr39,(n.B). 



(65) 



Focusing on the two crucial (growing or decreasing) 
asymptotic terms in the odd-parity metric, we can now 
compare the definition of b(, namely 

^ oc <^ + bi i^-^j r-' \ YtrniO, ip), (66) 

to the stationary limit of the general gravito-magnetic 
fields in a local A-frame (see Eqs. (2.19) and (4.16) of 
Ref. [7]), 



Ba=Ba+ B+ 



(67) 



Inserting now 5*^ — giHj^, contracting Ba with 
n°- , and recalling that one has n"9a = dr and 
dL{r-^) = i-Yi^i - l)!!r-(^+i), one finds 



Hl 



(68) 



so that 

^3 



L- 



(69) 



The comparison with Eq. (66) finally yields 
4(£ + 2) 



— {2i-l)\\Gcji^biy-^ J 



(70) 



Gcj, 



1 bp 



fGM\ 

i{i+2) {2i-i)\\ y^) 

l(£-h2) (2£- 1)!! 



2e+i 



(71) 
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where we used the notation = c^^'^^hi of Eq. (61). As 
announced, we see that, modulo a numerical coefficient, 
the odd-parity analog of the i?-scaled Love number ki 
is the combination c^^"*~^5f, Eq. (61). In the odd-parity, 
quadrupolar case {£ ~ 2) the above link reads 



G<J2 = ^i2R'' 



1 , /gmY 



(72) 



while the explicit expression of j2 



C'b2 reads 



32 



5 (2c (12(y + l)c4 + 2{y - 3)c^ + 2{y - 3)c^ + 3{y - 3)c - 3?; + 9) + 3(2c - l)(y - 3) log(l - 2c)) 



(73) 



r 



VI. COMPUTATION OF THE SHAPE LOVE 
NUMBER he 

We have indicated above (following a recent study of 
the gravitational polarizability of black holes [24]) how 
one can generalize to a relativistic context the shape tidal 
constant (the "first Love number" ) introduced in New- 
tonian theory. Let us first point out that, though in 
general there is no direct connection between hi and kt, 
there is a simple relation between them in the case where 
the deformed object is a ball of (barotropic) perfect fluid, 
treated in Newtonian gravity. Indeed, in the Newtonian 
theory of tidal deformation (see [42, 43]), we have the 
result that an external "disturbing" tidal potential 



U' 



disturb 



t 



(74) 



deforms the constant-pressure (and constant-density) 
level surfaces (p — p{a)) of a fluid star into 



r{a) 



fi{a)Pi{cose) 



(75) 



Here /^(a) = {Sr/r)^ satisfies the Clairaut equation (in 
which p(a) indicates the mean density within < r < a) 



6p(«) 
p{a) 



(76) 



and is related to the disturbing tidal coefRcient ci via 

{2i + l)c£A^+i 



G{£ + m)M ' 



(77) 



where A is the surface value of a (i.e. A ~ i?, the undis- 
turbed radius) and where rjf = [0'fi{O')/f£{o-)]a=A denotes 
the surface logarithmic derivative of fi{a). The latter 
quantity is related to the "second" Love number ki via 



2{i + 'ni) ■ 



(78) 



On the other hand, using the definition of hi, i.e. 



5R 



MA) 



'gm/r' 



we find 



2e + i 
J+m' 



(79) 



(80) 



By eliminating r/i between Eq. (78) and Eq. (80), we 
finally get a simple relation between hi and ki, namely 



hi = l + 2ki. 



(81) 



For instance, a Newtonian 7 = 2 polytrope has a den- 
sity profile p{r) cx sinx/x, where x = irr/R, from which 
one deduces, using either the Clairaut equation, or the 
Newtonian limit of the {(. = 2) even-parity master equa- 
tion (27), namely [6] 

H" + -H'+( AnGp^ -1)h = 0, (82) 
r \ dp 7'^ J 

that H cx x~^^'^ J5/2ix). [Note that this result is mis- 
printed as a;+^/^ 75/2(2;) in [6]]. This leads to [6] 



^2^(7 = 2) = -^ + -^=^ 0.25991 



1 



15 



and therefore 



15 

2 



27r2 



1.51982. 



(83) 



(84) 



When generalizing the definition of hi to the relativistic 
context it seems that wc lose the existence of a func- 
tional relation between hi and ki. Let us explicate the 
meaning and implementation of the relativistic version 
of hi, Eq. (16). First, we define 5R/R as the fractional 
deformation of a sphere, embedded in an auxiliary 3- 
dimensional Euclidean manifold, such that the inner ge- 
ometry of this deformed, embedded sphere is equal to the 
inner geometry (induced by the ambient curved space- 
time) of the real, tidally-deformed neutron star surface 
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(considered at fixed coordinate time). In general, one 
would need to consider the Gaussian curvatures of both 
surfaces to express their identity (as done in the black 
hole case [24]). Here, things are simpler because we are 
using the Regge- Wheeler gauge. In that gauge, it is eas- 
ily seen that the inner metric of the surface r = r{9, ip) 
of an (even-parity) tidally deformed star is 



Rn 



6r 

1 + ( 2 K 

r 



(85) 



where r{9,(p) = Ra{l + 6r/r) is the radial coordinate 
location of the star's surface. [Here we absorbed the 
Yimi^T^) factors in K and 5r/r.] Because this inner 
metric is conformal to the usual sphere 5^ of unit ra- 



dius dn'^ = d0^ + sin^ 



it is easily checked that it 



would be the inner geometry of a flat-embedded sphere 
(linearly) deformed by 5R/ R, with 



5R Sr 1 
-R^--2^- 



(86) 



We can further compute the value of the coordinate de- 
formation 5r/r by using the fact that the logarithmic en- 
thalpy h = dp/ {e+p) must vanish on the star surface. 
Since h{r) = ho{r) + Sh{r), where ho{r) is the enthalpy of 
the background undcformed star, and where 5h = —H/2, 
we then find 



Sr ~ 



1 ^H 



2 \ h 



r=R 



1 f e+p 

~7~ 



H 



(87) 



r=R 



where the prime denotes dj dr. Finally, we have the "flat- 
equivalent shape deformation" : 



8R 



V 



rp 



H - K 



r=R 



As we said above, the metric variable K can be expressed 
as a linear combination [31] of H and H' , namely K = 
aiH' + a2-ff. The coefficients ai and a2 evaluated on 
the unperturbed surface r = R of the star read 



2cR 



ai 



a2 



{£-l)(£ + 2)' 



(89) 



£{£ + 1) 



4c^ 



1 - 2c 



- 2(1 - 2c) 



(90) 



In addition, using the TOV equations, we also have on 
the star surface 



rp 



e+p I -2c 



(91) 



By replacing (on the surface) also RH' — yH{R), we 
finally obtain 



6R 



1 



H{R) 



1 - 2c 



2cy 



Q!2 



(92) 



2"'"' { c {£-!){£ + 2) 
At this stage, we have obtained an expression of the form 
6R/R = H{R)f{c,y). To proceed further and compute 
hi, it remains to obtain the value of JJ^"^*"''''. Following 
Ref. [24], and the spirit of the Newtonian definition of 
hi, we define ~2Uf^^*^^''^{R) as being the analytic con- 
tinuation at radius r = i? of the leading asymptotically 
growing piece in H, i.e. the part of ffsrowmg _ QpP^^ 
which grows as r^. In other words, we define it by the 
Newtonian-looking formula 



disturb 



(i?) 



-ap 



(93) 



We can then compute Uf^^*'™^{R) in terms of the full 
value of H on the surface H{R) — apPe2{x) + aQQi2{x) 
(with, we recall, x = R/AI — 1 = 1/c — 1) by separating 
two "correcting factors" out oi H{R), namely 



'-H{R)^ c'Pi2ix) 



1 + CLi- 



{12 



{x) 



' Pn{x) 



J /-disturb 



(94) 

Putting all the pieces together, and inserting our general 
result (39) for at — aq/ap, we obtain the following final 
result for he: 



hi = c'+' Pn{x) 



I -2c 



1 



{£-l){£ + 2) 



2cvi + £{£ +l) + 



4c^ 



1 - 2c 



2(1 - 2c 



1 - 



dx ^ogPnjx) - cyi 
dx logQnix) - eye 



= l/c-l 

(95) 



VII. RESULTS FOR THE EVEN-PARITY TIDAL the dependence of these quantities on the compactness 
COEFFICIENT /i^. 



Having explained how to compute the various tidal 
constants fig (or ki), ug (or jt) a-nd hg, let us discuss 
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c = GM / (cpi?), for various kinds of EOS. As we already 
mentioned, we shall consider a sample of EOS. 



A. Polytropic Equations of State 

First, we consider two kinds of relativistic polytropes: 
the "energy-polytrope" , or e-polytrope, such that p = 
Ate''', where e is the total energy density, and the "rest- 
mass-polytrope" , or ^-polytrope, with p = np? and 
e = pi-\-pl (7 — 1), where /j, = nrub is the baryon rest-mass 
density. For these polytropes, we shall focus on the adi- 
abatic index 7 = 2, which is known to give a rather good 
representation of the overall characteristics of neutron 
stars. We shall also briefly explore what happens when 7 
takes values larger (or smaller) than 2. In particular, we 
shall discuss, in the next subsection, the limit 7 00, 
which leads to an incompressible model, with uniform en- 
ergy density e. Let us note that the compactness of the 
7 = 2 e-polytrope models ranges between (for a for- 
mally vanishingly small central pressure) and 0.265 for 
the maximum mass model, while the compactness of the 
/Lt-polytrope ranges between and 0.2145. Note that the 
limit of a vanishing compactness c = GM/ (cqR) 
formally corresponds to the Newtonian limit. Augment- 
ing 7 allows one to reach higher compactnesses, and 
thereby to better explore the effects of general relativis- 
tic strong-field gravity. In particular, the incompressible 
limit, 7 — )■ 00, yields a range of compactnesses which ex- 
tends up to Cmax = 4/9 = 0.4444 . . . , quite close to the 
"black-hole compactness" c^^ = 1/2 = 0.5. We note 
that a theorem guarantees that 4/9 is the highest possi- 
ble compactness of a general relativistic perfect fluid ball 
(see e.g. [44, 45]). 

Figure 1 exhibits the dependence of the dimension- 
less, even-parity. Love number ki on the compactness 
c of 7 = 2 polytropes for three values of the multipole 
order: £ = 2, 3 and 4. The results for the /i-polytrope 
(solid lines) are compared with those for the e-polytrope 
(dashed line). The limiting values of kg for c 
(which numerically means c ~ 10~^) do agree, as they 
should, with the known Newtonian results [46, 47]. In 
particular, fe^ (7 = 2) = 0.25991 (as mentioned above) 
fcf (7 = 2) = 0.10645, and fcf (7 = 2) = 0.06024. 

The most striking structure of Fig. 1 is the very strong 
decrease of ki with increasing compactness. For typical 
neutron star compactness, say c ~ 0.15, the general rela- 
tivistic value of ke is about 4 times smaller than its New- 
tonian estimate k^ . This might have an important nega- 
tive impact on the measurability of neutron star charac- 
teristics through gravitational- wave observations. Leav- 
ing this issue to a future investigation [48], let us focus 
here on a deeper understanding of the c-sensitivity of ki. 

The main origin of the strong decrease of kt, when c 
increases, is the universal presence of an overall factor 
(1 — 2c) ^ in kg. As discussed above, the presence of this 
term is linked to the no-hair property of black holes, i.e. 
the fact that the star-rooted contribution cx aqQn in 




FIG. 1: Polytropic EOS; gravito-electric Love numbers kt 
(or apsidal constants) for £ = 2, 3, 4 versus compactness c = 
M/R. We use two different polytropic EOS's, either of the 
rest-mass type (p = np?; solid lines) or of the energy type 
(p = Ke^; dashed lines). For both EOS's we use 7 = 2. Note 
that the maximum compactness allowed by the e-polytrope is 
larger than that for the /i-polytrope. 



the metric variable H around a tidally-deformed star be- 
comes singular, in the black hole limit R 2GM/c^. In 
addition, there are other c-depcndent effects that tend 
to decrease the value of ki. This is illustrated in Fig. 2 
which plots the ratios kg = kg /[k^ {1 — 20)^] for the 7 = 2, 
/i-polytropic EOS. In the case of k2, Fig. 2 shows that the 
"normalized" Love number k2 is, to a good approxima- 
tion, a linearly decreasing function of c, k2 (c) ~ 1 — /3c, 
with a slope /3 ~ 3. In other words, the c-dcpendcnce of 
k2 (for 7 = 2) is approximately describable as 



k2{c) 



k^{l-2cf{l 



■(3c) 



(7 = 2), 



(96) 



with /? ~ 3. To get a more accurate representation, one 
must include more terms in the c-expansion of the "nor- 
malized" ^2, or more generally kg, say 



kg = fcf (1 



2c) 2^a„c 

n=0 



(97) 



Such a nonlinear fit yields an extremely accurate repre- 
sentation of the c-dependence of kg{c). The performance 
of such fits is illustrated in Fig. 2 (dashed lines) and the 
best fit values of the coefficients a^, < n < 4 (fitted 
over the full range < c < c„iax) are listed in Table I. 
Note that, if we were to trust these fits beyond the range 
< c < Cinax(7) where kg{c) is defined, they would pre- 
dict that ki{c) vanishes for a value cl slightly smaller 



than 



„BH 



1/2, and would become negative before van- 



ishing again (now quadratically) at c = c^^. The critical 
value c^ is around 1/3, and approximately independent 
of ^. 
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FIG. 2: Normalized Love numbers versus c for a 7 = 2 /x- FIG. 3: Incompressible EOS: gravito-electric Love numbers 

polytrope (points); and performance of the fitting with the ki (or apsidal constants) for ^ = 2, 3, 4 versus compactness 

template given by Eq. (97) (dashed lines). The coefficients of c = M / R. 
the fit for each I are listed in Table I. 



TABLE I: Fitting coefficients for ki as defined in Eq. (97) for 
a 7 = 2 /i-polytrope, up to £ = 4. 



I 


2 


3 


4 


ai 


0.9991 


0.9997 


0.9998 


ai 


-2.9287 


-5.0933 


-7.1938 


ai 


-1.1373 


7.2008 


18.9509 


ai 


14.0013 


1.0826 


-21.8488 


ai 


-50.9711 


-18.7750 


4.9031 



B. Incompressible Equation of State 

To further explore what happens for large compact- 
nesses, we have studied in detail the limit 7 ^ 00, i.e. 
the incompressible EOS, e = const. Let us recall that in 
this case the TOV equations can be solved analytically 
giving 



P 



VI - 2c - Vl - 2cr2 



Vl - 2cr2 - 3V1 - 2c 



(98) 



Here r denotes the dimensionless ratio r^^^^/R, so that 
< r < 1. Note that the central values (r ~ 0) of the 
pressure are pc = e{y/l — 2c— 1)/ (l — 3-y/l — 2c) . so that 
Pc 00 when c 4/9, which shows that c^ax = 4/9 is 



2c-l)/(l-3Vr^ 
4/9, which shows that c„ 
the maximum compactness reachable by an incompress- 
ible star. 

Let us now discuss the computation of ki for an in- 
compressible star. The limit e = const creates a tech- 
nical problem in the use of the master equation (27). 
Indeed, the coefficient Cq of H, Eq. (29), contains a con- 
tribution cx (e + p)de/dp which formally vanishes in the 



incompressible limit e = const. However, as e(r) is an 
inverted step function which vanishes outside the star, 
e{r) = Co (1 — 0(r — 1)), this term actually contributes a 
term oc d{r — 1) which must crucially be taken into ac- 
count in the computation of fc^. We can then proceed 
as follows. First, one numerically integrates the incom- 
pressible limit of Eq. (27) in the open interval < r < 1 
representing the interior of the star. The output of this 
integration is the value of the logarithmic derivative of 
H at r ~ 1~, say ym = y{R~). Second, one corrects this 
value into the value yout = y{R~^) just outside the star. 
To compute the correction let us evaluate the "strength" 
of the delta-function singularity in the only singular piece 
of Co, namely Cg'"^ = '^nGe^ie + p)de/ dp. Using the 
TOV equations giving the radial derivative of p, we can 



write (reverting to general units, with r ~ 

^sing ^ 47rG'r2 de 

m{r) + AnGr^p dr 



rPhys 



(99) 



In the incompressible limit de/dr = —eoS{r — R), m(r) 
(47rG/3)eoiT!^ and p{R) = 0, so that 



Cr' = +^Sir-R) 



(100) 



Then the effect of the singular term in Eq. (27), or, more 
clearly, in the corresponding Riccati equation for y{r) = 
rH'/H, 

ry' + y{y-l)+rCry + r^Co = 0, (101) 

is easily found to introduce a step function singularity in 
y{r) with strength y^^'^^{r) = ~39{r — R). This shows 
that the correct value ?/out = viR^) to be used in evalu- 
ating kg is (independently of the value of £) 



,,out 



yi 



3. 



(102) 
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As a check on this result, we can consider the Newtonian 
hmit of Eq. (27). In this hmit, the exterior solutions 
Pi2{x) and Qi2{x) reduce to {r/MY and (A//r)^+^ re- 
spectively, so that one has 

^ (103) 



1 



y 



which generalizes the i — 2 result of [6] to an arbitrary 
I. Then, the incompressible limit, in the interior, of the 
Newtonian limit of Eq. (27) reads 



H" + -H' 



H = 0, 



(104) 



which coincides with the exterior, Newtonian equa- 
tion for H, with general solution H = ap{r/MY -f 
aQ(M/r)^+^. Regularity at the origin selects the ap 
term, so that (r) = £ = vTi^)- Then, Eq. (102) 
determines 



so that 



out N 

Vl 



(incomp) 



4(£-l)- 



(106) 



This result agrees with the known result for a limiting 
7 = oo (n = 0) polytrope [46]. 

Figure 3 shows our numerical results for the ki Love 
numbers of a relativistic incompressible star, as function 
of c, in the range < c < Cmax = 4/9. We exhibit the 
three first multipolar order £ = 2,3 and 4. Note that the 
c — !■ values of k£{c) agree with the Newtonian limit, 
Eq. (106). The phenomenon of the "quenching" of fcf 
as c increases is even more striking, in this incompress- 
ible case, than in the 7 = 2 case considered in Fig. 1 
above. Note, in particular, the very small values reached 
by ki for the maximum compactness c = 4/9. To under- 
stand better the "quenching" of ki by strong-field effects, 
we have analytically studied the incompressible model in 
the limit of maximum compactness c 4/9. This limit is 
singular (because Pc — > 00), but is amenable to a full an- 
alytical treatment of ki. This analytical study is thereby 
a useful strong-field analog of the analytical study of the 
7 = 2 model in the weak- field (Newtonian) limit. Let us 
sketch how one can analytically solve the 7 = oo,c^4/9 
model. First, by introducing the variable 



Equation (27) becomes 



+ (4a; 
-f- {2a 



3 ' x^-Ax 



£(£ + 1) 



(107) 



dH 

dx 

-1)H = 0. (108) 



We found two exact, analytical solutions of this equation, 
which are both of the form (1 — x)°'{l -\- x)^ with some 
rational exponents a, (3. More precisely, either 



1 



a 



1 



(109) 



a 



e + 2 



(110) 



2 ' ''2 

Regularity at the origin ^ selects the first solution, namely 

H{x) = (1 - x)(^-i)/2(l + x)-(^+i)/^ (111) 

As a result, the interior value of the logarithmic deriva- 
tive j/int of H has the simple form 



y{x) = 1. 

X 



(112) 



Adding the effect of the (5- function at the surface ( r = 1, 
i.e. X = 1/3), Eq. (102) finally leads to 



yint(l/3)-3 = 3£-4. 



(113) 



(105) Inserting this result in our general result (49) for ki gives 



an analytical expression for /j,(™'^°"^P) (c^ax) • For instance, 
one finds the following value for £ = 2 

4096 



.(incomp)/ N 
'^2 V^niaxy 



10935(308- 81 log 3) 



0.0017103. 



(114) 

Note the striking quenching of fc2, by nearly 3 orders of 
magnitude, from the c ^ value /c^ = 0.75 to this result 
for Cmax =4/9. 



C. Realistic Equations of State 

Before discussing our results for the other tidal coeffi- 
cients {bf , <Ji, hf ), let us end this section devoted to the kf 
(and fj.£) Love numbers by briefly considering the kg, for 
the dominant quadrupolar order £ = 2, predicted by the 
two reahstic (tabulated) EOS FPS and SLy^. We recall 
that we have chosen them here because they have been 
used in some recent numerical relativity simulations of 
coalescing neutron star binaries [28, 29]. The maximum 
compactness of the two "realistic" EOS that we retained 
are c^^l = 0.2856 and c^Z, = 0.303. 

The corresponding results for kg are shown in the left 
panel of Fig. 4 for 2 < £ < 4. The right panel focuses 
only on the results for the £ = 2 case, that are plotted 
together with several illustrative /x-polytropes , namely 
7 = 1.8, 2 and 2.3. The 7 = 1.8 polytrope illustrates the 



* Note, however, that because of the singular behavior of p{r) at 
the origin, the "regular" solution H{r) is less regular than usual: 
H{r) oc r^~^ instead of r*. 

^ Since these EOS are given through tables, to use them in a nu- 
merical context it is necessary to interpolate between the tabu- 
lated values. As in [27] we use simple linear interpolation (instead 
of third-order Hermite or spline ones) to avoid the introduction 
of spurious oscillations in the speed of sound. See Ref. [27] for 
further details. 
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FIG. 4: The gravito-electric Love numbers ke_ (or apsidal constants) for 1 — 2,3,4 versus compactness c — M/R for the two 
tabulated "reaUstic" equation of state FPS and SLy (left panel). Right panel: comparison between k2 from various relativistic 
^-polytropes (with different 7) and the FPS and SLy realistic EOS's. 



reason why the "realistic" EOS lead to a decrease of kg 
as c ^ 0. This is Bacause of the fact that the "local" 
adiabatic index 



1 



dlogp dlogp 
d log e d log /i 



(115) 



of these EOS varies with the density (or pressure). As 
shown, e.g., in the bottom panel of Fig.l of Ref. [27] for ~ 
nuclear densities it stays in the range 2 < F < 2.3, while, 
for lower densities (around neutron drip) it drops to low 
values F < 1, before rising again towards F ~ 4/3 for low 
densities. Let us recall in this respect that Newtonian 
polytropes have a finite radius only for 7 > 1.2 (n < 5). 
When 7 — > 1.2, a Newtonian poly trope has finite mass, 
but its radius R tends to infinity. As k2 uses a scaling 
of ^2 by a power of R, this causes fc^ = ^2(0 = 0) to 
tend to zero as 7 ^ 1.2 (see [46]). Anyway, the decrease 
of fc2°^'"''"^(c) as c 0, linked to the small value of F 
for low (central) densities and pressures is a mathemat- 
ical property which is physically irrelevant for our main 
concern, namely the tidal properties of neutron stars. In- 
deed, neutron stars have a minimal mass determined by 
setting the mean value of F equal to the critical value 
~ 4/3 for radial stability against collapse [49]. More- 
over, we are mainly interested in neutron star masses 
~ 1.4Mq. Such neutron stars are expected to have radii 
varying at most in the range 10km ^ R ^ 15km, corre- 
sponding to compactnesses 0.13 ^ c < 0.2. To be on the 
safe side, we shall consider the interval 0.12 < c < 0.22, 
which is indicated by vertical lines in the right panel of 
Fig. 4. Focusing our attention on this interval, we can 
draw the following conclusions from the inspection of the 
right panel of Fig. 4: 



bracket the ^2(0) sequence predicted by the two 
realistic EOS's. 

2. Actually, the two realistic EOS's that we retained 
here, FPS and SLy, lead to rather close predictions 
for fc2(c). 

3. In the range 0.12 < c < 0.22 the two, "reahstic" 
fc2 (c) can be approximately represented by the fol- 
lowing linear fit 



(116) 



with A ~ 0.165 and B ~ 0.515. 



VIII. RESULTS FOR THE ODD-PARITY TIDAL 
COEFFICIENTS ae 

As explained above, the odd-parity tidal coefficients 
is obtained by solving the master equation (31). Let us 
start by noticing that the formal Newtonian limit of this 
master equation is simply 



(117) 



Indeed, all the matter-dependent contributions to 
Eq. (31) arc, fractionally, of order Gm{r)/{cQr) ^ 
47rGr^e/cg or AnGr'^p/cQ, and vanish in the Newtonian 
limit c —^ 0. In this limit, Eq. (117) does not contain any 
effect of the star, and, in particular, is the same in the 
interior or in the exterior of the star. This shows that 
the origin-regular solution of Eq. (117) is, everywhere, of 
the form 



1. The /^-polytropes 7 = 2 and 7 ~ 2.3 approximately 



■0 (r) oc r 



(118) 
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FIG. 5; The j2, odd-parity, Love number for the 7 = 2, /i- 
polytrope and for the incompressible EOS. 



and does not contain any "decreasing" , Q-type contribu- 
tion bqtpQ oc r^^ This proves that the odd-parity Love 
number be ~ bq/bp vanishes in the Newtonian Umit. 
More precisely, as the first post-Newtonian corrections 
to Eq. (117) are fractionally of order c = GM/{cqR), 
we can then easily see that the i?-normalized odd-parity 
Love number jg = c^^~^^bg will be of order c as c ^ (in 
agreement with the results of [23] concerning the 1 = 2 
case). 

As we have pointed out above that bi also contains a 
factor 1 — 2c, we conclude that j£ vanishes both when 
c ^ and c — > 1/2, and should qualitatively be of the 
type 



31 = c''+'be ~ Bec{l - 2c). 



(119) 



This approximate result suggests that a IPN-accurate 
calculation of the solution of the master equation (31) 
should give us access to the coefficient Bi, and thereby, 
in view of Eq. (119), to a global understanding of the 
c-dcpcndcncc of je. In particular, one expects from (119) 
that = |c^^+^6f| will attain a maximum value some- 
where around c = 1/4, with the value 



\B,\ 



(120) 



Therefore, a IPN computation of the coefficient B gives 
an indication of the maximum strength of the odd-parity 
Love number. We have analytically computed the co- 
efficient Bi (defined by jg = Bgc + C(c^)) by solving 
Eq. (31) by perturbation theory ip = '0o + V'l- Here 
ipo = r''+^ is the solution of the c ^ Hmit, Eq. (117), 
of Eq. (31), and tpi is the first-order effect of the matter 
terms in Eq. (31). Actually, we found convenient to get 
tp = Ipo + ipi by using Lagrange's method of variation of 
constants: ipi''') = Ci(r)r^"'""'^ + C2{r)r~^ , with Ci(r) — > 1 



and C2(r) ^ as r ^ 0. Skipping technical details, we 
found that the logarithmic derivative 



y 



rip' 



(^+1) 



1 



£+1 ci(r) 



(121) 



takes the following value at the star surface 

e + 2 1 



2/(i?) = (£ + l) 1 + 



X / drr 
Jo 



£ + 1 i?2£+l 

2e 



2{£-2)— +A'Kr\e-p) 
r 



(122) 



The integral term in this result represents the IPN cor- 
rection Bacause of the presence of matter. On the other 
hand, the small-c limit of the general result (61) yields 
for £ = 2 (chosen for simplicity) 



32 



ILzl 

y + 2 



y- 



(123) 



where we used the fact y = 3 + 0{c). Combining this 
result with Eq. (122) above yields 



J2 



/ drA-wGie - p)r^ . 
Jo 



(124) 



One can analytically compute this integral in the case of 
a Newtonian polytropc with 7 = 2. Recalling that in this 
case we have e — p ~ pcQ ~ PcCq sinx/x, with x = Trr/R, 
the result is 



j2=S2C + 0(c2) (7 = 2) 



with 



Bo 



4 / 20 120 

5 \ 7r2 ^4 



-0.164395, 



(125) 



(126) 



in agreement with [23] (modulo normalization issues that 
we did not check). 

In view of the reasoning above, we then expect that 
j2 — B2c[\ — 2c) will be negative, will vanish at c = 
and (formally) at c = 1/2, and will reach a minimum 
value (j2)"=2 - B2/S = -0.02055 around c = 1/4. For 
completeness, let us also mention that the case of an 
incompressible neutron star leads to 



B. 



(incomp) 



= -0.342857 (127) 



and (j2)™eomp - -0.042857. We have qualitatively 
and semiquantitativcly confirmed these results on the c- 
dependence of the odd-parity Love number by numeri- 
cally integrating Eq. (31). We display in Fig. 5 the result- 
ing odd-parity quadrupolar Love number j2 , versus c, for 
both a 7 = 2 /i-polytropc (solid line) and an incompress- 
ible EOS (dashed line). We have numerically checked 
that the slope at the origin of the c-axis is indeed B2 as 
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FIG. 6; The j2, odd-parity, Love number for FPS and SLy FIG. 7: Shape Love numbers hi versus c for the two 7 = 2 



EOS. 



polytropic EOS: the /^-polytrope (sohd hnes) and the e- 
polytrope (dashed hnes). 



analytically determined above. In both cases (though it 
is more evident in the incompressible case, where higher 
values of c are allowed) j2 has a negative minimum before 
rising again towards zero. The numerically determined 
minimum of ^2 is min(j2) ~ —0.0216 (reached around 
c ~ 0.21) for 7 = 2 and min(j2) — —0.0634 (reached 
around c ~ 0.315) for the incompressible EOS. 

From the conceptual point of view, these results on the 
odd-parity Love number (and, via Eq. (72), on the cor- 
responding magnetic-like tidal coefficient ai) are inter- 
esting counterparts of the even-parity results discussed 
above. They have points in common (their vanishing in 
the formal limit c 1/2), and they also strongly differ 
in other aspects: vanishes when c 0, while /i^ has 
a well-known Newtonian limit; and ui is proportional to 
the first power of 1 — 2c, while /i^ oc (1 — 2c)^. More- 
over, cTf, as naturally defined, is negative, while /i^ is 
positive. As, with our DSX-like normalization, the in- 
teraction energies associated to both types of couplings 
arc proportional (modulo positive numerical constants) 
to —MlGl = —iiiG\ and —SlHl = — cr^iJl respec- 
tively, this sign difference can be interpreted as being 
linked to the well-known, Lorentz-signature related, fact 
that current-current interactions have always the oppo- 
site sign to charge-charge, or mass-mass, interactions. 
Concerning the formal vanishing of cr^ as c ^ c^^ = 1/2, 
the same remarks we made above for the even-parity case 
apply here. This fact is essentially, given the no-hair 
properties of black holes, a consistency check on the def- 
inition of fJi as measuring a violation of the "effacing 
principle" . As said above, although it suggests that the 
correct value of ui for black holes might be zero, it is far 
from proving such a statement which has a meaning only 
within a more complex nonlinear context. 

From the practical point of view, an interesting output 
of the investigation of (T2 is that its numerical value hap- 
pens to be quite small. Indeed, for a 7 = 2 /i-polytrope 



we have 

I _ I max , 

= — hsT'''^ ~4 X 10"^ (128) 

We shall discuss in another work the precise dynamical 
meaning of this small number (e.g. for the dynamics of 
binary neutron stars) , but the appearance of such a small 
number clearly means that it will be an enormous chal- 
lenge to measure it via gravitational- wave observations. 

For completeness, we conclude this section by showing, 
in Fig. 6, the behavior of j2 also for the two different real- 
istic EOS, FPS and SLy, that we have introduced above. 
Similarly to the case of fc2 (see Eq. (116) above), in the 
range 0.12 < c < 0.22 (i.e., between the two dashed ver- 
tical lines) the two "realistic" j2 (c) can be approximately 
represented by the following linear fit 

_^.(FPS;SLy)^^_^^^ (129) 

with ^ ~ 1 X 10"" and B ~ 0.1411. 



IX. RESULTS FOR THE "SHAPE" LOVE 
NUMBERS hi 

Equation (95) gave the final expression for the (even- 
parity) "shape" Love number hg. This expression con- 
tains several terms that are singular as c c^^ = 1/2: 
(i) the long curly bracket contains an explicit term oc 
1/(1 — 2c); (ii) the logarithmic derivative of (3^2 be- 
haves, when X —* \ (given that Qi2{x) ^ [x - l)-i), 
as dx\ogQn{x) ~ —{x — 1)"^; and (iii) the logarith- 
mic derivative of Pf2 behaves, when a; ^ 1 (given that 
Pg2{x) ~ [x — 1)^^, sec below), as dx\-OgPn{x) ~ 
+ {x — 1)~^. The latter behaviors mean that, in the 
limit R/{2M) 1, the last factor in Eq. (95) tends 
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FIG. 8: Shape Love numbers hi versus c for the incompress- 
ible modeL 



to 1 + 1 = 2. As for the (1 — 2c) singularity in the 
curly bracket, it is compensated by the linear vanishing 
of Paix) as X 1^ indeed 



Pt2{x) ^{x" -1) 



1 



1-2. 

c 



(130) 



Finally, the formal "black-hole limit", c -> c"" = 1/2 of 
the "shape" Love number hi{c) is finite and nonzero. Ac- 
tually, we found that this limit agrees with the results of 
a recent direct investigation of the "gravitational polar- 
izability" of a black hole [24] . The general result for hf^ 
can be found in the latter reference. Let us only mention 
here the values of the first two "shape" Love numbers: 



lim /12(c) = 



BH 



lim /13(c) = /i?" = — . 

3V ^ 3 20 



(131) 
(132) 



Figure 7 shows the results of inserting the numerically 
determined value of ye{R) into the expression (95) of hi. 
We give the results for the first three multipolar orders, 
£ = 2,3, 4, and for the two 7 = 2 polytropes (^-polytrope 
and e-polytrope). We have also investigated the results 
for the incompressible EOS (7 — *■ 00). They are shown 
in Fig. 8. This information is completed in Fig. 9, where 
we investigate the effect of changing the EOS on the c- 
behavior of the leading, quadrupolar, shape Love number 
/i2. In all cases, we sec that, somewhat similarly to the 
/c2 case, the strong self-gravity of a neutron star tends 
to "quench" the value of /i2. For instance, as we dis- 
cussed above, the Newtonian limit of a 7 = 2 polytrope 
yields h^{-/ = 2) = IS/tt^ ~ 1.52. As we see in Fig. 9 
this value is reduced below 1, i.e. by more than 33%, 
for typical neutron star compactnesses. When exploring 
stronger self-gravity effects, notably for the incompress- 
ible model, one gets an even more drastic quenching of 



FIG. 9; Influence of the EOS on /12(c): comparison between 
the incompressible case and ^-polytropes with two values of 
the polytropic index: 7 = 2 and 7 = 3. The red circle on the 
right of the plot indicates the formal c — > 1/2 result. 



/i2, by an order of magnitude, from the Newtonian value 

J A'' (incomp) 



value /if" 



= 2.5 down to a value near the "black hole" 
= 1/4 = 0.25. From the theoretical point of 
view, it is nice to see this continuity, as the compactness 
increases, between the neutron-star case and the black- 
hole case. Note that neither the no-hair property of black 
holes, nor the related "effacing principle" , are relevant to 
the present result. What is relevant is that the inner ge- 
ometry of the horizon of a black hole is well defined and 
that a black hole is an elastic object, like a neutron star. 

Similarly to what we did for discussed above, it 
is also convenient and useful to represent hi as a c- 
expansion of the form 



n=0 



(133) 



where h^ is the Newtonian value (obtained from 
through Eq. (81)) and the coefficients 6^ are obtained 
from a fit. As an example. Table II lists these coeffi- 
cients for a 7 = 2, yu-polytrope up to £ = 4 (i.e., they are 
obtained by fitting the solid lines in Fig. 7). 

For completeness, we conclude this section by dis- 
cussing the hi results for the two different realistic EOS, 
FPS and SLy, that we have introduced above. In Fig. 10 
we display the hi Love numbers (for £ ~ 2,3,4) versus 
compactness c. The fact that hi ^ 1 when /c^ — > (be- 
cause of the small value of the local adiabatic index F for 
low central densities and pressures) is understood via the 
Newtonian link (81). 
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FIG. 10: Shape Love numbers he versus c for the two tabu- 
lated realistic Equation of State FPS and SLy. 



TABLE II: Fitting coefficients for he as defined in Eq. (133) 
for a 7 = 2 /i-polytrope, up to ^ = 4. 



£2 3 4 



&o 0.9999 0.9999 0.9999 

bi -3.6764 -4.3700 -5.2361 

bi 4.5678 6.9775 10.4578 

bi -0.0192 -4.1964 -9.6026 

bi -5.8466 -1.028 3.0415 



X. CONCLUSIONS 

We have studied the various tidal responses of neu- 
tron stars to external tidal fields. We have considered 
both electric-type (even-parity) and magnetic-type (odd- 
parity) external tidal fields. As indicated by Damour, 
SofFel and Xu [7] some time ago, one can correspond- 
ingly introduce two types of linear response coefficients: 
an electric- type tidal coefficient Gfii = [length] mea- 
suring the i^^ mass multipole GMl induced in a star 
by an external ^'^-order (electric) tidal field Gl, and a 
magnetic-type tidal coefficient Gai = [length]^^+^ mea- 
suring the i*"^ spin multipole GSl induced in a star by 
an external .^^'^-order "magnetic" tidal field Hl. Divid- 
ing Gfig and Gcr^ by the {2i + l)-th power of the star's 
radius R leads to dimensionless numbers of the type in- 
troduced by Love long ago in the Newtonian theory of 
tides. In addition, one can define a third^'^ dimensionless 



Love number (for any £) , measuring the distortion of the 
shape of the surface of a star by external tidal fields. 

We have studied, both analytically and numerically, 
these various tidal response coefficients, thereby gener- 
alizing a recent investigation of Flanagan and Hinderer. 
The main results of our study are: 

1. A detailed study of the strong quenching of the 
electric-type tidal coefficients fie (or its dimension- 
less version kg ~ Gfie/R^^^^) as the "compactness 
c = GM/{cqR) of the neutron star increases. This 
quenching was studied both for polytropic EOS (of 
two different types, see Fig. 1), for the incompress- 
ible EOS (where the quenching is particularly dra- 
matic, see Fig. 3) and for two "realistic" (tabu- 
lated) EOS (see Fig. 4). 

2. Part (though not all) of this quenching mechanism 
can be related to the no-hair property of black 
holes. The latter property ensures that some of 
the tidal response coefficients of neutron stars must 
vanish in the formal limit where c —> c^^ = 1/2. 
At face value, this suggests that the "correct" value 
of the fj,e and ae tidal coefficients of black holes is 
simply zero. We, however, argued that this conclu- 
sion is premature, until a 5PN (5-loop) nonlinear 
analysis of the effective worldline action describing 
gravitationally interacting black holes is performed. 

3. We gave accurate nonlinear fitting formulas for 
the dependence of the tidal coefficients kg and he 
of a 7 = 2 /x-polytrope on the compactness (see 
Eqs. (97) and (133)). We also found that two "re- 
alistic" EOS give rather close values both for the 
electric and magnetic tidal coefficients of neutron 
stars. In particular, this suggests a possible, ap- 
proximately universal analytical representation of 
the leading, quadrupolar (electric) Love number 
for neutron stars of the expected compactnesses, 
0.12 < c < 0.22, namely 

fc2(c) ~ 0.165 -0.515c. (134) 

Even if this simple linear fit reproduces with only 
a few percent accuracy the c-dependcncc of the 
known realistic EOS, it might suffice to deduce from 
future gravitational-wave observations an accurate 
value of the neutron star compactness. Indeed, the 
dimensionless parameter which is crucially entering 
the gravitational-wave observations is the dimen- 
sionless ratio 

The strong c-dependence of /i2(c) coming from the 
power implies that even an approximate fit such 



One could also introduce magnetic-like "shape" Love numbers 

by considering other aspects of the geometry around the star surface. 
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as (134) might allow one to deduce from the mea- 
surement of /i2(c) a rather accurate (say to better 
than 1%) cstimate^^ of c. 

4. We surprisingly found that the magnetic-type Love 
numbers of neutron stars are negative, and quite 
small. We showed, by analytical arguments, that 
they can be approximately represented as oc Bc{l — 
2c) with a calculable coefficient B (that we com- 
puted in a few cases). 

5. Following a recent investigation of the gravita- 
tional polarizability of black holes [24], we stud- 
ied the "shape" Love numbers hi of neutron stars. 
Again the quantity hi{c) is found to be drastically 
quenched when c increases. However, in that case 
hg{c) does not tend to zero as c c^^ = 1/2. 
Rather we found that hf^{c) tends to the nonzero 
black-hole value hf^ [24] as c formally tends to 
c™ = 1/2. 

In future work, we will come back to the other issues 
mentioned in the Introduction, namely: 

1. the incorporation of tidal effects within the effective 
one body formalism, starting from the additional 
term in the effective action 

AS = +^^2 / dsf.^f (136) 



2. the study of the measurability of various tidal co- 
efficients within the signal seen by interferometric 
detectors of gravitational-waves. 

After the submission of this work, a related paper by 
Binnington and Poisson [50] appeared on the archives. 
Ref. [50] develops the theory of electric and magnetic 
Love numbers in a different gauge. Their results seem 
to be fully consistent with ours, but are less general: (i) 
their treatment is limited to e-polytropes, (ii) they did 
not consider the "shape" Love numbers, and (iii) they do 
not discuss the effective action terms associated to tidal 
effects. 
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